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Abstract 

There are now evidences that the cosmological constant A has a 
non-zero positive value. Alternative scenarios to a pure cosmologi- 
cal constant model are provided by quintessence, an effective negative 
pressure fluid permeating the universe. Recent results indicate that 
the energy density p and the pressure p of this fluid are constrained by 
—p < V ~ —0.6 p. Since p = —p is equivalent to the pure cosmological 
constant model, it is appropriate to analyze this particular, but im- 
portant, case further. We study the linear theory of perturbations in a 
Friedmann-Rober-tson- Walker universe with a cosmological constant. 
We obtain the equations for the evolution of the perturbations in the 
fully relativistic case, for which we analyze the single- fluid and the 
two-fluid cases. We obtain solutions to these equations in appropriate 
limits. We also study the Newtonian approximation. We find that 
for a positive cosmological constant universe (i) the perturbations will 
grow slower in the relativistic regime for a two-fluid composed of dark 
matter and radiation, and (ii) in the Newtonian regime, the pertur- 
bations stop growing. 
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1 Introduction 



Ever since Einstein first introduced it, the cosmological constant A has had 
a history full of ups and downs. The interest in it has been reawakened 
by novel results, which now lead us to believe that A is an essential fac- 
tor in the present dynamical evolution of the universe. A result supporting 
the existence of a cosmological constant is due to the recent balloon-borne 
measurements of the cosmic background radiation ( ||de Bcrnardis et al. 2000 



Hanany et al. 20U0] ]). These measurements lead to a non-zero cosmological 
constant and to an Q ~ 1 universe. Another important result comes from the 
study of high redshift supernovae ( fKiess et al 1998| , [Ferlmutter et al. 1999| ), 
which leads us to believe that the expansion of the universe is accelerat- 
ing. In combination, both studies yield for the present energy densities 
of matter and the cosmological constant roughly the values Qm — 0.3, 
Q\ ~ 0.7. An alternative to the cosmological constant scenario is pro- 
vided by quintessence, a scalar field furnishing an effective negative pres- 
sure fluid ( Peebles 1984| , Peebles fc Ratra 1988|| ). Recent results indicate 
that the energy density p and the pressure p of the fluid are constrained by 
-p < V ~ -0-6 p IIEfstathiou 20001 . Since p = — p is equivalent to the pure 
cosmological constant model, it is appropriate to study this case further, in 
some detail, which we shall do in this paper. 

A long standing issue yet to be fully solved, in which the cosmological 
constant A (or variants such as quintessence) can play a role, is how the 
multitude of structures we observe in the sky has formed from a homoge- 
neous background. In harmony with the novel findings mentioned above, 
some works have also shown that the existence of a positive cosmological 
constant in a cold dark matter scenario yields results that fit better with 
observations, namely, the power in the perturbation spectrum as a func- 
tion of wavelength (or mass) ( ||Kofman et al. 1993| , |Krauss fc Turner 1995 



Ustrikcr fc Stcinhardt 1995|j , and the cluster abundance as a function of 



redshift | |Bertschinger 199"B| |Bahcall et al. 19TJ9] ]). There are now some good 
reviews on the effect of a non-zero cosmological constant on the universe 
( ICarroll et al. 1992j |Cohn 1998; , |Sahni fc Starobinsky 20001 , |Carroll 2000] ). 

One subject that has not been discussed in detail for a universe with a 
positive cosmological constant, and which is a basic step in the understand- 
ing of structure formation, is the study of the equations and solutions of 
linear perturbations from the homogeneous background. This has been done 



2 



exhaustively for models considering mainly dark matter, both warm and 
cold, normal matter and radiation (see, e.g., [Peebles 1980| , Peebles 199 3 



Padmanabhan 1993| ], and references therein), but mostly leaving out the ex- 



istence of a possible cosmological constant. Some exceptions are the papers 
by [[Kofman fc Starobinskii 1985|| where the authors study anisotropies in the 
microwavebackground in a positive A universe, ||Silveira fc Waga 1994|| where 
a perturbation analysis and power spectrum calculation in a universe with a 
time-dependent cosmological constant is performed, [[Abramo fc Finelli 2001 



and [[Riazuelo fc Uzan 2000|1 where one finds a perturbation analysis in mod- 
els of quintessence, and [[Proty 2000|| in which an extrapolation from a A = 
case to a A > case is made. 

An excellent discussion on the theory of linear perturbations and their 
connection to structure formation is given in the book of Padmanabhan 
( [[Padmanabhan 1993 |). Our work presented here is based mainly on chapter 
4 of Padmanabhan's book. It is our aim to extend the analysis and study of 
the linear perturbation theory by including the effects of a positive cosmo- 
logical constant. A Newtonian analysis has been presented in ( ||Cohn 1998|| ), 
which is valid for perturbations with wavelengths much smaller than the 
horizon size and for non-relativistic fluids. We enlarge the analysis to the 
relativistic regime, valid for wavelengths larger than the horizon size, and for 
relativistic fluids. As it is known, the existence of a positive A affects the rate 
of expansion of the universe, but not the the self-gravity of the perturbations. 
This causes a suppression on the growth of the perturbations. Since obser- 
vations favour a flat universe we set Q = 1 in our study. We then find that, 
indeed, for a positive A-Universe (i) the perturbations will grow slower in the 
relativistic regime for a two-fluid composed of dark matter and radiation (the 
growing mode grows with the power 1/4 instead of 1 which characterizes the 
A = case), and (ii) in the Newtonian regime, the perturbations stop growing 

T9931 |Cohn 1998|| ). 



Kofman et al 



in a A-dominated universe (see also 

Our paper is delineated as follows. In section 2 the perturbed equations 
for a single relativistic fluid are found and solved for a radiation dominated 
phase, for a matter dominated phase and for a A dominated phase. In section 
3 the perturbed equations for two relativistic fluids are studied and applied for 
the coupled system of radiation and dark matter. In section 4 the perturbed 
equations in the Newtonian approximation are obtained and solved. These 
are valid for length scales smaller than the horizon size. In section 5 we 
conclude. 
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2 The relativistic case: single fluid 



2.1 Equations 

When studying perturbations in general relativity there is always the prob- 
lem of fixing a gauge, or coordinate system. A discussion of the pertur- 
bation equations in different gauge conditions is given by Hwang and Noh 
( ||Hwang & N oh 19971 , see also |[Bardeen 19801 ), while Unruh ( ||Unruh 1998)1 ) 



gives a discussion on the problems of gauge fixing. We will fix our gauge as 
the comoving gauge, with the four-velocity u a given by u a = (1,0,0,0), to 
derive the perturbation equations for a nonzero A universe. The derivation of 
the equations for A = is given by Padmanabhan ([ [Fadmanabhan 1993| |); we 



will here follow his procedure closely, modifying the equations to introduce 
the cosmological constant. 

We use the overdot to denote the derivative u a D a . The projection tensor 
onto surfaces orthogonal to comoving world lines, the comoving hypersur- 
faces, is h a p = g a p — u a up. The natural derivative on these comoving hy- 
persurfaces is h^Dp, and the Laplacian is given by V 2 = —bP a Dph° n D 1 . The 
metric is the usual Friedmann- Robert son- Walker metric for a flat universe, 
given by 

ds 2 = dt 2 - a{t) 2 [dr 2 + r 2 {d9 2 + sin 2 Odcj) 2 )} . (1) 

The energy momentum tensor is T a/3 = (p + p)u a u^ — pg 01 ^, where p is the 
energy density and p the pressure of the fluid. We start with the relativistic 
continuity and Euler equations, derived from D a T^ = 0: 

p=-3H(p+p), (2) 
«.= ^. (3) 

P+P 

The Einstein field equations, G a p — Ag a p = 87iGT a p, yield the Friedmann 
equation 

Note that it is sometimes interesting to consider the cosmological constant 
term, Ag a @, as being an energy- momentum tensor for a special fluid. Here 
we do not follow this procedure, since we find more useful to see it as a 
cosmological constant term. 
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Expanding ii a } using the definition for overdot given above, we have 

D a ii a = (D a u p ) (D p u a ) + u l3 Dp{D a u a ) + u p R a(3 u a . (5) 

In our chosen frame, the first two terms in the left hand side are purely 
spatial. Using 

u a R a(5 u p = R 00 = 4vrG(p + 3p) - A , (6) 

and, splitting the tensor into a diagonal term (which in the absence of 
perturbations will be related to the Hubble's constant H by Hubble's Law) 
and non-diagonal terms (which are perturbative terms whose square is of 
second order) we get, up to first order 

D a ii a = 3H + 3H 2 + 4vrG(p + 3p) — A. (7) 

Using now Euler's equation (|]), and again up to first order, we can relate 
the left hand side of (0) to V 2 p by 

D a u a = -^y. (8) 

p + p 

Combining equations (0) and (|8|), we have 

H = -^-^ {p + 3p) + ^ ^L. (9) 
3 3 p + p 

Following Padmanabhan ( fPadmanabhan 1993|| ), we need to change our time 
variable to have a valid time label for the comoving hypersurfaces. Taking t 
to be the valid ordering label, we have: 

± = i-JE- (io) 

dt p + p 

Taking now the perturbed quantities of the form A = Af, + 5 A, where Af, is 
the background, non-perturbed quantity and 5 A is a small perturbation, we 
can rewrite the continuity equation (|2|) up to first order as 

¥ = p b + 5p+^^5p (11) 
dr p b + p b 
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where the overdot is used to represent the derivatives with respect to t. 
Taking the linearized right-hand side of @, and equating the zeroth order 
terms, we get 

p b = -3H b {p b + p b ), (12) 
and, for the first order terms, 

Sp = -35H{p b + p b ) - 3H b 5p . (13) 

For H, we have 

^- = H b + 5H b + -^5p. (14) 
dr p b + p b 

Linearizing equation (H) and equating to equation (|i~4|), we have for the zeroth 

order terms, 

H b = -H 2 b -^(p b + 3p b ) + ±. (15) 

Using now the zeroth order expression fll5|) and Friedmann's equation, (|4|), 
the first order term is 

S H = -2H > SH- 4 ^S P -i^, (16) 
3 3 p b + p b 

where we used the definition 5p = c 2 6p, c s being the speed of sound. We 
can then rewrite equation ( fL3|) in first order in 5H, and differentiate it with 
respect to t to get: 



SH 



3 (1 + w) 



6 - 3H b {2w -c 2 s )5+ ^H%(2cl -w + w 2 )5- -Aw(l + w)5 



(17) 

where w = p/p and 5 = 5p/ p, and we have used 

w = -3H b (l + w)(c 2 s - w) } (18) 

H b = J-Hl{l + w)+ l -K{l + w). (19) 

Using (|I~3"D, ( |I~6D and (p^7|), we finally arrive at: 

5 + H b (2 -Qw + 3c 2 s ) 6 - \H 2 {\ - Qc 2 - 3w 2 + 8w)5 

mi+w)(l-3w)A5 = -(^) 2 5, (20) 

where we introduced the Fourier transform such that V 2 5 = —(k/a) 2 5. The 
Newtonian limit, is obtained by taking w ~ and c 2 ~ in the left hand 
side and using Friedmann's equation M) (see also section [§). 
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2.2 Solutions 

We start by changing the independent variable in equation ( p0[ ) from t to 



a: 



J 2 5 



d8 



<r-— r + Aa— 

da z da 



k 2 r 2 



where 



B = -2(1 - 6c s 2 - 3w 2 + 8w) + 2(1 + «;)(! - 3^)^ 



(21) 

(22) 
(23) 



Radiation dominated phase: 

We take radiation to be the dominant component, and study its perturba- 
tions Sr. We take (k 2 c 2 )/(H 2 a 2 ) <C 1, which in this case corresponds to con- 
sidering perturbations larger than the Hubble radius. Taking c 2 = w = 1/3, 
equation (^l|) becomes 



2 d 2 5n t n A d§R 



da 2 



+ 2 



3H 2 da 



-25 



R 



0. 



(24) 



We now rewrite Friedmann's equation as 

8%G 



^(Pr + Pdm) + ^ • 



(25) 



Ignoring the matter term in (^5), we can write H as H = Hq(£Ir a\ja +0/ 



where the quantities VL are the usual quotient of the density to the critical 
density, given by £Ir = ffffPi? and = j^. Using this in equation (|24]) , 
we get the solutions 



S Rg oc a 2 2 F 1 



rl i 7 -- 



Q\ a 4 



5Rd 



4' '4' n R a* 
JIao 4 + £IrOq 



(26) 



(27) 



where d and g stand for the decaying and growing modes and 2-F1 [a>, b; c; z] 
is the hypergeometric function which can be given in terms of a power series 
Y^h=o( a )k(b)k/ (c)k z h /k\. As f^a 4 <C QrOq, these solutions reduce to those 
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obtained in the A = case: 5r 9 oc a 2 , 5nd oc 1/a, as would be expect from 
the fact that, in the radiation dominated phase, the cosmological constant 
does not have much influence in the dynamics of expansion, given by equation 

(§• 

To consider perturbations smaller than the Hubble radius in the radiation 
dominated phase, we take (k 2 c 2 )/(H 2 a 2 ) 3> 1, and ignore the B term in 
(|2T1). As we are considering the radiation dominated phase, the cosmological 
constant does not have much influence in the dynamics of evolution, and 
thus the solution in this case is the same as the one for the A = case (see 
Padmanabhan 1993|[ ). 



Matter dominated phase: 

We now study the perturbations in dark matter, 5 dm in the matter domi- 
nated phase (to study matter perturbations in the radiation dominated phase, 
we need to consider the two fluid situation, developed in the next section). 
We have c 2 = w = 0, and, using this (but preserving c s in the {k 2 c 2 )/(H 2 a 2 ) 
term to account for pressure support), equation ( pip becomes 

2 d 2 5 3 / A \ dS f k 2 c 2 s 3 3 A \ n 
a da~ 2 + 2 ( X + MP^a + (h4 ~ 2 + 2 MP^ = ° ' (2§) 

We take (k 2 c 2 )/(H 2 a 2 ) <C 1, which here corresponds to considering wave- 
lengths larger than the Jeans wavelength, Aj oc Cg/p 1 ^ 2 , which is the char- 
acteristic scale at which pressure support becomes an important factor. For 
matter domination, we approximate (p5|) by taking only the matter term. 
We can then ignore the A terms in (p8]), and get the usual solutions (see, 
e.g., fPadmanabhan 1993|| ), given by 5dm = Aa + Ba^l 2 . 



Lambda dominated phase: 

If we take into account the cosmological constant, as well as matter, we 
can write H as H 2 = Hq(Qm «o/« 3 +^a), where VLm = tjpPdm- Substituting 







this into equation (^8f) and once again ignoring the term {k 2 c 2 )/(H 2 a 2 ), we 
will get the solutions: 



■3 



6 ^ d = A—^- (29) 
6 DMg = Bx 2 F 1 [-,1,^, -x 3 ] (30) 
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where 2F1 is the hypergeometric function defined earlier, and x is a mod- 
ified scale factor defined by 1 = a/a\ = xa eq /a\, where x is the usual 
scale factor x = a/a eq , with a eq the scale factor at the time of matter ra- 
diation equivalence and a a = a (fi M /^A) 1//3 the scale factor of the universe 
when the energy density of dark matter and of the cosmological constant is 
the same. These solutions have also been obtained by Silveira and Waga 
[Silvcira fc Waga 1994 ] for a study of perturbations in a universe with a time 



dependent cosmological constant. The behaviour for zero cosmological con- 
stant is obtained by taking the limit 1 < 1, and corresponds to the solutions 
obtained for the matter dominated phase. Solution (^) is shown in Fig. 1. 



Figure 1. Amplitude of the dark matter perturbations as a function of the 
modified scale factor x; the dashed line represents the evolution of the per- 
turbation for A = 0, the full line in the presence of A. 

One sees that for f < 1 the perturbations follow the behaviour for A = 
0, 5dm oc x, while for x 1, when the cosmological constant dominates 
dynamically, the growth of perturbations is suppressed. 

We finally consider the case where the cosmological term dominates in 
equation (f|). We can then write H 2 ~ A/3 and equation ( p8| ) becomes simply 



2 d 2 5 DM d5 DM „ / 01 x. 
a— —- + 3— — = 0, (31) 
aa z da 

whose solution is 

S DM = A + Ba- 2 . (32) 

Comparing this to the solution for the matter dominated phase, given by 
5 = Aa + Ba~ 3 / 2 , we can see that the presence of the cosmological constant 
causes the perturbations to stop growing. 
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The growing mode throughout the different phases is then given by: 

5r = a 2 (radiation d.p., mode bigger than the H.r.) 

* 5dm = a (matter d.p., mode bigger than J.l.) 

5dm = constant (A d.p., mode bigger than J.l.) 

where "d.p." stands for dominated phase, "H.r." for Hubble radius, and 
"J.l." Jeans length. The one fluid relativistic description is good enough to 
study perturbations in the dynamical dominant component. However, if one 
wants to understand the behaviour of the other fluid components one needs 
a two-fluid description. This we shall do now for the case of radiation and 
dark matter. 



3 The relativistic case: two fluids 

3.1 Equations 

In order to study simultaneously the perturbations in radiation and dark 
matter, we shall first set up the equations for the relativistic, multi-fluid 
case, and then specify for a 2-fluid with those components. We again follow 
( IPadmanabhan 1993|| ), generalizing the procedure to include a non-zero cos- 



mo logical constant. 

We consider the case of a perfect fluid made up of N perfect fluids, each 
with an energy momentum tensor given by Tff = (px + Pn)u n Un — pNg a/3 , 
where, as before, we work in the comoving gauge, u N = (1,0,0,0). Each of 
these fluids will then obey equations (Q) and (|3|), where the quantities will 
now have indexes N. Following the same procedure as in the single fluid 
case, we get the following equations: 

p N = -3H N (p N + p N ) , (33) 
±D a u% = H N + H 2 N + ^(p + p) - A , (34) 



W 2 p N , 3(H-H 



D a u N = ■ 1 ■ p N . (35) 

Pn +Pn Pn+ Pn 

In equation (|3~5"|), the last term in the right hand side is new when compared 
to equation (H). This comes from the fact that, in the equation for each fluid 
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equivalent to (|3]), we no longer have h£, which is the correct projection tensor 
onto the comoving hypersurfaces, but a term h N ^ = — UN a u N . This new 



term reduces to zero in the case of a single fluid. We also have for the total 
fluid quantities: 



P 



Y.PN 



p 



N 



N 



^PN+PN u 

N 



P + P 



and 



H-H 



N 



5H-5H 



N ■ 



(36) 

(37) 
(0) and 

(35), using the change in time variable described in the last section, given by 
fllOD , and equation (|37|) , to get the following perturbation equations: 

5p N = -3(p N + p N ) - 3H5p N - 3H(5p N - 9 N 5p) (38) 

Pn 



We can now follow the procedure of linearizing equations 



where 



4vrG x 1 V 2 p N 
—^—OP - o : 1- 







3 pn +Pn Pn + Pn 
Pn + Pn 



Mj 



M 



N 



P+P 



SH N 
(39) 

(40) 



The quantities wn and defined before for the fluid N. For the 

total fluid these are given by: 



P 

w = - 
P 



cl = V - 



p 



(41) 



where p and p are defined by (|36|). It is important to notice that, unlike what 
happened with a single fluid, we have here c\ ^ 5p/5p. The final equation 
for the perturbations is then given by 



5 N + (2 + 3a 

rl 



2 



6wm)H5 



n 



3H — c s N 



1 



w 



M P 



■-sM 



-sat 



M 



+ 3 



-3H 2 + A)(c sN - w N )(l + w) + 5H z (c sN - w N ) + 2Hc sN c : 



sN 



>N 



2 l + ^ ( 2 

+ 1 + w {s 



w 



2 

, 2 PM ( C&M 



~ W M 



M 



CI — W 



r 2 
S N 



>M 



11 



6H2 1 h -3# 2 + A))(l + w N ) ^2 —c^ N 8 M 



Vl+w 2 > ^ p 

3H +Wn ^ E*£.qjj ( w _ Wm ^ + 2c sM ) - 4vrG(l + w N )5 
1 + w y p 



-SAT 



V 2 5iv = (42) 



We are interested in studying the two fluid case, which in this case will 
be radiation and dark matter. Following ( [padmanabhan 19931 ), we now 
introduce the adiabatic 5 and isocurvature S perturbations related to the 
individual fluid perturbations by: 

5 = Plh + M2 ; S = ^ (43) 

Pl + p2 l + Wi 1+W 2 

Using these definitions, the two perturbation equations become: 
5 + H(2 - Qw + 3c 2 )<5 - §# 2 (1 - 6c s 2 + 8w - 3w 2 )<5 

+§(l + w)(l-3w)A5 = -§(c s 2 5 + w?) (44) 

and 



where 



S + H(2 - 3« 2 )S = ^((c s 2 - c s 2 )-^- - u 2 s) , (45) 
er \ 1+w ) 

(46) 

pif(l + U7) 
2 = Pl(l + Wi)c s | + p 2 (l + W 2 )C S 2 

P(l+W) ' 1 ' 

Equation ( fH|) reduces to that of the single fluid case ( p0|) if we take r/ = 
and a single fluid component. 



3.2 Solutions 

We discuss here the case of a two component system, with radiation and 
dark matter. We will use as variable, instead of a, x = a/a eq , with a eq being 
the scale factor at the time of matter-radiation equivalence. With this new 
variable we can write: 

Pdm = 1_. £* = ±_. ^ = ^L±PH = l j{x + lh (48) 

Peq Peq Peq Peq 
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Taking wdm = = c s 2 DM and wr = 1/3 = c sR , we have 



w 



PDM + PR 



PDM + PR 



Pdm + Pr 3(1 + x) ' s pdm + Pr 9 x + 4/3 
Using these definitions, we may write H as 



H 2 (x) = J-;(x + wi + 4(1 



2x 4V ~ ' eq ' 3 
We will characterize the perturbations by the parameter u, with 

k 2 2x 2 , 



H 2 a 2 



-LU 



X 



(49) 



(50) 



(51) 



with to related to the wavelength by 2-kuj = [dji(t eq ) / X(a eq )], where dn is the 
Hubble radius. Perturbations for which u > 1 will enter the Hubble radius 
in the radiation dominated phase, while for u < 1 they will enter it in the 
matter dominated phase. 



Using the new variable x, and definitions (^8[) to 
x(d/dx), we can rewrite equations Q44]) and (f45|) as 



L> 2 5 + 



5 x 



g _ 1 + I A f Q 
2 1+x x+4/3 ^2 3H 2 1 u 



1 

1+x 



D5 



II), and defining D 



1 x 2 i 3x 

2 (1+x) 2 " t " 4 



I 9 x' 2 _ 3x^_ _ o , n x+4/3 A W 
r 4i+4/3 1+x- Z ZX (x+1) 2 3H 2 1 



3 (x+l)(x+4/3) 



(xS - (JC + 1)5 



2 x 2 

3 x+1 



2x 4 



A 

i ; 3// 2 



(52) 



D 2 S + 



2(l+x) 



I 1 OTX 

"""3I+4/3 

= 1 ^ 
3 x+4/3 



x+4/3 
^2 



+ 



1 



A 





2x 2 N 


) a" 




x+l y 


1 3/F 



1+x / 3H 2 



5 



2x 2 + ( 2£1 



2x 2 ) 3H2 



5 



(53) 



The perturbations 5 and 5* can be related to the perturbations in the radia- 
tion and dark matter components by: 



3 x + 1 

A R ~ x + 4/3* 



x 



x + 4/3 



S. 



(54) 
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x + 1 4 i 

x + 4/3 3x + 4/3 

With these relations we can now write the equations in Ar and Adm corre- 
sponding to equations (|52| ) and fl5"3"|): 



£ 2 + =: 



x 



3 x + 4/3 A 
2 x + 1 + 2 aT+1 3# 2 



6 + ? 



2 Z 



-(— 

1 V3F 2 



A 



DA/ 



3 x + 4/3 



D + 2 



a; + 4/3 



A 



R 



l + x-2x 4 



A 



9 x 2 (x + l)(x + 4/3)3H 2 



u 2 A } 



(56) 



f)2 i /l i,3 x+4/3 A \ f) i ( 2 A! , 4 1 / 

\ 2 x+1 2 x+1 3.ff 2 / V 3 :E+1 3 z+4/3 Vx+4/3 



1 2 1 2i 2 i 4 x 3 \, ,2 A 



"•"l x 2 3 x(x+4/3) x+1 "r" 3 (x+l)(x+4/3) r° 3H 2 



A 



3 x / i A I x 7-) 

2 x+1 I 1 3H 7 j x+4/3 



Ada/ • (57) 



We are interested in studying the solutions of these equations in the 
two limiting cases u 2 ^> 1 (perturbations entering the Hubble radius in the 
radiation dominated phase) and w 2 C 1 (perturbations entergin the Hubble 
radius in the matter dominated phase). We can further divide our study to 
the cases where x C 1 and x>l. As we have discussed earlier, the effects 
of the presence of a cosmological constant will only be important for x ^> 1, 
well into the matter dominated phase, as the cosmological constant starts 
to assume a greater importance in the dynamics of evolution. Therefore, we 
will here restrict ourselves mainly to these cases; we will, however, discuss 
a situation where i C 1 in order to certify that the solution obtained is 
approximately the same as that of the case A = 0. 

uj 2 ^> 1: perturbations entering the Hubble radius in the radiation dominated 
phase 

The case uo 2 3> 1, for which the perturbations enter the Hubble radius in 
the radiation dominated phase, has three distinct cases: (i) when the mode 
is bigger than the Hubble radius, ux <C 1, the solution is very close to that 
of the A = case, given by Ar = Ada/ = x 2 , so we will not discuss it 
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further here; (ii) when ujx ^> 1, but i < 1, the mode has entered the Hubble 
radius, but the universe is still in the radiation dominated phase; this case 
also shows very little deviation from the zero cosmological constant case, but 
we will consider it here to show this; (iii) finally, there is the case ujx ^> 1, and 
x ^> 1, when the universe is in the matter dominated phase; this case shows 
significant deviation from the A = case. For zero cosmological constant, 
( ||Naycri fc Padmanabhan 1998 ]) give an interesting approach to this case, 



based on a Newtonian theory with a modified continuity equation (see also 



Lima et al. 1997 



We now discuss case (ii). Approximating equation ( |57| ) for ujx ^> 1, 

x <C 1, we get 

2 uj 2 

D 2 A R - DA R + -uj 2 x 2 A r + —n A A R = . (58) 

Using equation (0), we can write ^a as: 

A 2x 4 

Qa = mp = QA ^ x + i + n Acci (2x*-x-i) (59) 

where the index eq refers to the time of matter radiation equivalence. This 
expression can, in this case, be approximated to Q A — ^Aeq2x 4 , where we 
took Q Aeq <C 1. As ^Acq is constant, equation (^) can be solved to give 

Ar = Aexp (±iux) ; v 2 = (- + 2Q Aeq )u; 2 > 1 . (60) 

This represents a rapid oscillation at frequency z/, which, due to the fact that 
^Acq <C 1, is very close to the frequency of oscillation of the solution for the 
A = case (see fPadmanabhan 19931 ). The solution for Adm, A dm = hix, 



is also approximately that obtained for zero A. 

We finally consider case (iii), x ^> 1. The equations are: 

t) 2 A DM + -(1 + 3n A )DA DM + -(ft A - 1)A DM = —uJ 2 n A A R , (61) 



+ 3n A )DA R + ^uj 2 x(1 - n A )A R = ~(1 - Q A ) A DM - DA DM . 

(62) 
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For a value of x sufficiently low, the perturbations will still follow the be- 
haviour they had in the A = long as the cosmological con- 
stant does not produce significant changes in the dynamics of expansion. 
From ( [[Padmanabhan 1993|| ), we see that the growing mode of A DM is Ax, 
while the corresponding mode of A R is (3A)/(4w 2 ). For ^a close to 1, the 
term in the first derivative in the left hand side of ([H]) is then close to 
2 Ax, while the term in the right hand side is 4v4/3; taking x ^> 1, we 
can then ignore the term in Ar. Substituting Q\ for A/3H 2 , and taking 
H 2 = Hq{VLm a o/ a eq l/^ 3 + ^a)j where the quantities Qm and Q\ now refer 
to their present values, we obtain for the solutions to the equation in Adm- 



x j/n 



^ K —2 > ( 63 ) 

tt M a 
-3' ' 6 I^a a o 



d g (xx 2 F 1 ,-— ^-x* . (64) 



These correspond to the solutions obtained previously to the equations for 
the perturbations in matter in the single fluid relativistic case, (|29|) and (|30|) . 
We now take the limiting case £7a — 1. In this limit, equations ( |61~D and (|S^) 
become: 

D 2 A DM + 2DA DM = , (65) 
D 2 A R + DA R + DA DM = . (66) 
The solution of ( |55D will then be 

A dm = Sdm = A + Bx~ 2 , (67) 
and, taking the dominant mode in (E7]), the solution of ( |56D is 



Ar = ~6 r = C + Dx- 1 . (68) 

As expected, the presence of a cosmological constant causes a halt in the 
growth of perturbations, with the perturbations in dark matter being greater 
than the radiation ones. 

We can summarize the evolution of the perturbations as: 

Sdm = x 2 ; 5r = 4/3x 2 before entering the H.r. 



'DM 

6 DM = lux 
5dm 



oscillates radiation d.p., after entering the H.r. 
8r tend to a const, for A domination 

after entering the H.r. 
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uj 2 <C 1: perturbations entering the Hubble radius in the matter dominated 
phase 

We now treat the case u 2 <C 1. In this case, the modes enter the Hubble 
radius in the matter dominated phase. The solutions for x <C 1 are approx- 
imately the same as in the case A = 0, (see fPadmanabhan 1993||). We find 



the solutions for the case w 2 x 2 3> 1 and x ^> 1. As there is a coupling 
(through terms of order ux) between S and S in equations fl52"D and (|53|), 
we can take 5 = S here, and, with these approximations, equations (|52[) and 
)3T) become: 

&S+(-l + lA)DS = 0. (70) 



v 2 ' 23H 2 , 

In the A dominated phase we can approximate f2 A = A/(3H 2 ) as a constant, 
since if will be approximately constant. The dominant growing mode is then 
given by: 

1 



6g = S g = x n °, n 9 = -(-l-3fi A + V 25 - 18fi A + 9fi i); (71) 



and the decaying mode by 



1 



5 d = S d = x n «, n d = -(-l-3ft A - V25 - 18fi A + 9fii) . (72) 

For fi A — 1; these exponents become n g = and = —2, which correspond 
to the solution fl32|) . For fi A = 0.7, we have n g = 0.25 and n d = —1.8. 
In the limit A — > 0, these reduce to the appropriate solutions, n s = 1 and 
= —3/2. The dominant modes can then be summarized as follows: 



2 



5 = x 
5 = x 
5 = x 713 



S = \nx before entering the H.r. 

S = x matter d.p., after entering the H.r. 

S = x™ 9 limiting case for A d.p. after entering the H.r. 



4 Newtonian Approximation 
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4.1 Equations 



The equations in the Newtonian approximation are derived classically, 
so the introduction of a cosmological constant will not change them; it will 
only influence the solutions when we consider the dynamics of expansion. 
We will here give a very brief description of the procedure and present the 
equations we will be using. We start with the classic fluid equations of mass 
conservation and pressure support: 

P=% + (v • V)p = -p(V ■ v) (73) 

v = -V0 - p _1 Vp (74) 

We relate the velocity v to the Hubble constant H by Hubble's Law, v(i, x) = 
H(t)x, and consider quantities to be of the form A = A b + 5A. We substitute 
these in equations ([73]) and (|74]), and linearize them, taking only the terms 
up to first order. Using V 2 <50 = 4irG5p and the definitions 6p = c 2 5p and 
5 = Sp/p, we arrive at the following equation for the perturbations: 

5 + 2H b 5 - c 2 s a- 2 V 2 5 = ^Gp b b , (75) 

where the a~ 2 term in the Laplacian comes from considering Friedmann co- 
moving coordinates. Considering perturbations with wavelength A, and tak- 
ing the Fourier transform, we can write the V 2 term as —k 2 . Equation (fT5|) 
then takes the form: 

5 + 2H b 5 + ^5 = AnG Pb 5 . (76) 
a 2 

This discussion can be generalized for the multi-fluid case. The corre- 
sponding equation is then 

5 N + 2H b 5 N - c 2 N a- 2 V 2 5 N = 4vrG ]T p M S M (77) 

M 

where the sum on the right hand side is over all components. A smoothly 
distributed component (e.g. the cosmological constant term) does not con- 
tribute to this term. 
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4.2 Solutions 



We consider here the case where A ^> Xj. We can then ignore pressure 
support, that is, ignore the k 2 term in the left hand side of ( [76|) , and get the 
equation 

6 + 2H5 - AnGpd = (78) 

As shown by Heath ( [[Heath 1977|| ), an integral expression for the growing 
solution of equation (|78|) for any value of f2 M and fl^ can be obtained by 
noticing that H{t) is itself a decaying solution to this equation. The growing 
mode can then be given by 



5 oc H(z) 



00 1 + x 
HHx) i ' 



(79) 



with 



H(z) = H (l + z) 



Q M z + ^a(q" 



1 



1/2 



(80) 



and z = a /a — 1 is the redshift. 

An analytical solution for the A = case is widely known (see, for exam- 
ple, [ Padmanabhan 1 993|| ) , and an analytical expression for ( |79|) in terms of 
elliptic integrals is given by Eisenstein ( [[Eisenstein 1997]] ). A good approxi- 
mation to this integral, given by 



5 



a -\l M 



4 / 7 o 

i4 M i4 A 



1 



2 Qm 



(81) 



(which is normalized to give 5 = a for = 1 , — 
Carroll et al. ( HCarroll et al. 1992|1 ). Taking the limit Q M 
the integral ([75|) can be approximated by 



S oc 



1 



1 - n. 



0) was given by 
- with Q A ^ 1, 



(82) 



which is a constant. We expect the perturbations to stop growing when the 
cosmological constant becomes the dominant energy form, because the expan- 
sion becomes faster than the gravitational collapse of the perturbations. The 



expansion timescale is given roughly by t exp ~ (Gpdominant) 1 ^ 2 ~ (Gpa)^ 1 ^ 2 , 
and is smaller than the collapse timescale, t grav ~ (Gpm)' 1 ' 2 , as long as 
Pa > Pm\ thus, the rapid expansion of the background prevents the growth 
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of the perturbations by gravitational collapse. 

We now study the perturbations in the matter component in a universe 
with a cosmological constant, in the case A ^> Aj where we can ignore the 
pressure support. The perturbation equation is simply (|78|), with H given 
by (p5|). In order to use the Newtonian equations, we are considering here 
the case when the universe is well into the matter dominated phase, and as 
such we can ignore the radiation term in (|25|). Further, as we are not taking 
into account the radiation, we need not distinguish between dark matter and 
baryonic matter, as they do not have different interactions with the cosmo- 
logical constant term, unlike what happened with radiation. 

Using the variable x = a/a^, we get the equation: 

(l + x 3 )x 2 5" + hx + 2x^)5 ' -~5 = (83) 

The solutions to these equations are the same as the ones obtained earlier 
for the relativistic case, given by (|29|) and (pCf). 



5 Conclusion 

We have studied the equations for the gravitational collapse of density pertur- 
bations in an expanding universe of the Friedmann- Robertson- Walker type, 
whose dynamics is described by equation fl25|), taking into account the pres- 
ence of the cosmological constant A (or equivalently, considering the universe 
permeated by a quintessence fluid with an effective equation of state given 
by P = —p)- We have confirmed that the cosmological constant does not 
have much influence in the behaviour of perturbations until it becomes a 
significant factor in the expansion of the universe, see equation (|25|). Thus, 
as we have explicitly shown for the particular case of perturbations in ra- 
diation for the two-fluid case in (|60|), the growth of the perturbations for 
A ^ is not significantly altered in the radiation dominated phase, when 
compared to a A = universe. On the other hand, when the cosmological 
constant becomes important in the expansion, at later times, we have seen 
that significant differences arise in the behaviour of perturbations. This is 
caused by the fact that the expansion becomes faster than the gravitational 
collapse of the perturbations, thereby preventing any further growth. This 
behaviour can be seen, for example, in the exponent n g of the solution (|7l|) , 
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which for A = is n g = 1, while for = 0.7 we have n g — 1/4. A more 
detailed solution for the transition from the phase of matter domination to 
that of A-domination has been obtained, both in the fully relativistic case 
and in the Newtonian approximation (which we have shown to be a good 
approximation in this case). This solution is given by equation (|30|) , and is 
illustrated in Fig. 1. In this figure we can clearly see the transition from the 
behaviour in the matter domination phase (x <C 1), when the perturbations 
grow as a, to the phase of cosmo logical constant domination (x ^> 1), when 
the perturbations stop growing altogether. 
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